function FeatureExtra_YCbCr(AU_DIR, SP_DIR, T)
    folderpathAu=AU_DIR;
    folderpathSp=SP_DIR;

    dstnmeAu=dir([folderpathAu '/*.bmp']);
    FmAu=size(dstnmeAu,1);
    single = (2*T+1)*(2*T+1);
    
    %% AU
    FeatureYAu_x = zeros(FmAu,single);
    FeatureYAu_y = zeros(FmAu,single);
    FeatureYAu_d = zeros(FmAu,single);
    FeatureYAu_i = zeros(FmAu,single);
    tic;
    for k=1:FmAu 
        fileNameAu=[folderpathAu '/' dstnmeAu(k).name];
        fprintf(1,'File : %s',fileNameAu);
        YAu = imread(fileNameAu);
        AU_DCT = abs(DCTImage(YAu));
        clear YAu;
        FeatureYAu_x(k,:) = markov_x(AU_DCT,T);
        FeatureYAu_y(k,:) = markov_y(AU_DCT,T);
        FeatureYAu_d(k,:) = markov_d(AU_DCT,T);
        FeatureYAu_i(k,:) = markov_i(AU_DCT,T);
        clear AU_DCT;
        fprintf(1,'\tFINISHED!\n');
    end
    time1 = toc;
    save(['my_FeatureYAu_x' num2str(T) '.mat'],'FeatureYAu_x');
    save(['my_FeatureYAu_y' num2str(T) '.mat'],'FeatureYAu_y');
    save(['my_FeatureYAu_d' num2str(T) '.mat'],'FeatureYAu_d');
    save(['my_FeatureYAu_i' num2str(T) '.mat'],'FeatureYAu_i');
    clear FeatureYAu FeatureYAu_x FeatureYAu_y FeatureYAu_d FeatureYAu_i
    
    %% SP
    dstnmeSp=dir([folderpathSp '/*.bmp']);
    FmSp=size(dstnmeSp,1);
    FeatureYSp_x = zeros(FmSp,single);
    FeatureYSp_y = zeros(FmSp,single);
    FeatureYSp_d = zeros(FmSp,single);
    FeatureYSp_i = zeros(FmSp,single);
    tic;
    for k=1:FmSp
        fileNameSp=[folderpathSp '/' dstnmeSp(k).name];
        fprintf(1,'File : %s',fileNameSp);
        YSp=imread(fileNameSp);
        SP_DCT = abs(DCTImage(YSp));
        clear YSp;
        FeatureYSp_x(k,:) = markov_x(SP_DCT,T);
        FeatureYSp_y(k,:) = markov_y(SP_DCT,T);
        FeatureYSp_d(k,:) = markov_d(SP_DCT,T);
        FeatureYSp_i(k,:) = markov_i(SP_DCT,T);
        clear SP_DCT;
        fprintf(1,'\tFINISHED!\n');
    end
    time2 = toc;
    save(['my_FeatureYSp_x' num2str(T) '.mat'],'FeatureYSp_x');
    save(['my_FeatureYSp_y' num2str(T) '.mat'],'FeatureYSp_y');
    save(['my_FeatureYSp_d' num2str(T) '.mat'],'FeatureYSp_d');
    save(['my_FeatureYSp_i' num2str(T) '.mat'],'FeatureYSp_i');
    clear FeatureYSp FeatureYSp_x FeatureYSp_y FeatureYSp_d FeatureYSp_i;
    %% output process time
    disp(['Au process time:',num2str(time1)]);
    disp(['Sp process time:',num2str(time2)]);
end